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Abstract 

An explicit relativistic light-front quark model is presented which gives the momentum 
transfer dependent form factors of weak hadronic currents among heavy pseudoscalar and vector 
mesons in the whole accessible kinematic region 0 < < q^ax ■ numerical investigations 

of the B —> 0*11^1, B ^ plvi, D ^ and D ^ plvi semileptonic decays the equal time wave 

functions corresponding to the updated version of the ISGW model are adopted. Using the 
available experimental information on branching fractions BR{B —> D*liyi) and BR{B —> plvi) 
the CKM parameters Uft and Vuh were estimated: lUfel = 0.036±0.004, \Vub\ = 0.0033±0.0004. 
The model is further tested by comparison with experimental data, QCD sum rules and lattice 
calculations. 


PACS numbers: 13.20 He; 13.20 Fc; 12.39 Ki; 12.15 Hh 

Keywords: semileptonic decays; heavy mesons, relativistic quark models 


1 


1 Introduction 


The knowledge and understanding of the form factors of hadronic currents is of decisive importance for 
the determination of the quark mixing parameters. In heavy-to-heavy transitions the determination 
of these form factors can be obtained based on heavy-quark symmetry in a nearly model-independent 
way [Q. On the other hand, in the theoretical description of the heavy-to-light decays there is 
no symmetry one can apply to constrain the relevant hadronic form factors. Quark model (QM) 
calculations can be very helpful for heavy-to-heavy as well as for heavy-to-light transitions 
Although strict theoretical error limits cannot be given, they provide a vivid picture of what is going 
on and give numerous testable predictions for quite different processes. 

In refs, p, ^ the dependence of the form factors among various pseudoscalar mesons has been 
calculated in the whole accessible kinematic region 0 < using the light-front (LF) QM. 

The calculations have been performed in a reference frame where the momentum transfer is purely 
longitudinal. The form factors were determined by using matrix elements of the plus components of 
the currents which are “good” operators in the LF formalism |^. The purpose of this paper is to 
generalize the results of refs, p, ^ and extend them for the calculations of the form factors governing 
pseudoscalar-vector {PS — V) transitions to 1“ states. We will show that the latter can again be 
expressed in terms of the matrix elements of the plus components of the currents in a reference frame 
corresponding to the purely longitudinal momentum transfers, but the calculations require the use 
of a limiting procedure. 

The plan of the paper is as follows. In Sections 2-4 we describe our approach for the calculation of 
the weak meson form factors. In Section 5 the numerical results for the form factors and the exclusive 
decay rates of the B D{D*)ii'i, B — 7i{p)ii'£, D K{K*)ii'£ and D tt{p)Ip£ weak decays 
are presented and compared with experimental data as well as with the results of other approaches. 
Section 6 contains a brief summary. 


2 Kinematics 


We denote by Pi, P 2 and Mi, M 2 the 4-momentum and masses of the parent and daughter mesons, 
respectively. The meson states are denoted as |P > for a pseudoscalar state and |P, e > for a vector 
state, where £ is the polarization vector, satisfying e • P = 0. The 4-momentum transfer g is given 
by g = Pi — P 2 and the momentum fraction r is dehned as 




r = 


= 1 - 


( 1 ) 


q 

Pt - pf 

In the following we work in the rest frame of the parent meson. The ^-momentum P 2 of the 
hnal meson is in the plane 1-3, so that q_L 7 ^ 0. We denote the angle between P 2 and the 3-axis by 
a. Then, it can be easily verihed that 


M| 


ql 


q^ = {l-r){Mi-^)-^-^, 
r r 


( 2 ) 


ql = Ml{y‘^ - (3) 

where the “velocity transfer” y is dehned as y = ui ■ U 2 with ui and U 2 being the 4-velocities of the 
initial and hnal mesons. The relation between y and g^ is given by 
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^ 2 M 1 M 2 ■ ^ ^ 

The momentum fraction of Eq. (1) is invariant under the boosts along the 3— axis and the rotations 
around this axis but depends explicitly on the recoil direction. Solving Eq. (2) for r one obtains 

r{q^,a) = C{y+\Jy^ - Icosa), (5) 

where ( = At the point of zero recoil r does not depend on a, = C,. 

The hadronic form factors for semileptonic decays are dehned as the Lorentz-invariant functions 
arising in the covariant decomposition of matrix elements of the vector and axial currents. We dehne 
the form factors of the PiiQiq) —^ P 2 {Q 2 q) transitions between the ground state A-wave mesons in 
the usual way. The amplitude < P 2 \V^\Pi >=< P 2 \Q 21 ^lQl\Pl > can be expressed in terms of two 
form factors 


<P2\V,\P,>= [P^ 


\ , Mf-Ml 


q^ 


-q^^ Fi{q) + 


q^ 


Qi^Paiq'^), 


where P = Pi + P 2 . There is one form factor for the amplitude < P 2 , V^l-Pi > 

2i 


<P 2 ,e\V^\Pi>= 




M 1 + M 2 


and three independent form factors for the amplitude < P 2 ,e\A^\Pi >=< P 2 ,£\Q 2 liil 5 Qi\Pi > 


( 6 ) 

(7) 


< P2-,£\Afj_\Pi > — ({Ml + M2)e*^Ai{q^) — ^ ^^ (Pi + P2)^A2{q^) — 2M2-^q^A^{q^) 


e*q 


+2M2^-g^Ao(g ), 


where Ao( 0 ) = ^ 3 ( 0 ) and the form factor A^{q^) is given by the linear combination 


A,{q^) = 




2 Mo 


2 M 2 


Sometimes another set of form factors is introduced 

< P2,£|V;|Pi >= tg{q^)e^uaqe*^P^q^ 

< P2,£|A^|Pi >= - [/(g2)e*+a+(g2)(PP)P^ + a_(g2)(PP)g, 
The relationship between two sets of form factors is given by 

V{q^) = -{Mi + M2)g{q^), 


( 8 ) 


(9) 

( 10 ) 

( 11 ) 

( 12 ) 


Ao{q^) 


\^[f{q^) + (^1 - M|)a+(g2) + q^a_{q% 


Ai{q^) 


fjq^) 

Ml + M 2 ’ 


(13) 

(14) 
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A2{q^) = (Ml + M2)a+(g^). (15) 

In the massless lepton limit the form factors FQ{q‘^) and a-{q^) do not contribute to the decay 
rates because of leptonic current conservation. The differential decay width for a semileptonic decay 
of a PS particle to another PS particle and a massless lepton pair {eVe or is given by 


dr G|Mf 4 2/ 2 1'|3/2| jt? 

( 2 / - 1 ) \FM\. 

where Vq^q^ is the CKM matrix element. The differential transition rate to a vector particle is 


dV GlMi 


dy 487r^ 
where the helicity amplitudes are given by 


C\VQ,Q,\\y^ - + H^), 


ddn = 


Mf 


{y-0(1 + 0Aliy)-2{y^-l)C 


My) \ 

i+cr 


= Ml (^(1 + C)dli(i/) T 2(1/2 - . 

Near the kinematic point of zero recoil, // = 1, Eq. (17) simplihes to 0 

while at opposite end of the heavy meson spectra, namely, at — 0 one obtains 


( 16 ) 

( 17 ) 

( 18 ) 

( 19 ) 

( 20 ) 


f = l^l'^«.«.tC(i-CV|A„(0)P. (21) 

In the case of heavy-to-heavy transitions, in the limit in which the active quarks have inhnite 
mass, all the form factors are given in terms of a single function Oy)i the Isgur-Wise form factor. 
In the realistic case of hnite quark masses these relations are modihed: each form factor depends 
separately on the dynamics of the process. 


3 The matrix elements of the vector and axial currents 

We determine all of the form factors by taking matrix elements of the plus components of the weak 
currents. Applying the standard LF technique for the time-like momentum transfer P, Q to the 
matrix element of the vector current ( 527 ^< 5 i between two pseudoscalar mesons one obtain^ 

Jviq\r) =< P 2 IQ 27 +Q 1 IA >= j^Sd^k^X2ix\k'l)xiix,kl)-Iv, (22) 

0 

where ly is the contribution of the Dirac current and quark spin structure: 

ly = /il/i2Tr[P+M(p2, \2)'l""u(py Ai)Poo], (23) 

^The spectator quark carries the fraction x of the plus component of the meson momentum, while the heavy quark 
carries the fraction 1 — x. 
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with 


and 


/ii = [M^o - (mi 

/i2 = [M|o - (m 2 - m)Y^^ 


M2o^M2o(a;,A;i) 


m? + k\ m\ + k\ 

^ i ’ 

X 1 — X 

w? + k'l m^ + 

x' 1 — x' 


(24) 

(25) 

(26) 
(27) 


In these eqnations mi and m 2 are the masses of active quarks, m is the mass of the quark-spectator, 
x' = ^ and k'j_ = k_L + x'q_L. In our kinematics 


k']^ = ki — x'q±, k '2 = k 2 - 


(28) 


The spin wave functions Rjj^ can be found in refs. H, We choose a phenomenological wave 
function Xi{^i ^1) by starting with the equal-time wave function Wi{k‘^) normalized according to 


fOO 

/ dk k'^ w‘l{k‘^) = 1, 

Jo 

provided that the fraction x is replaced by the longitudinal momentum k^^'^ dehned as 

1 ' 


k!'^^ = (X 


Mio -|- 


2 2 
m\ — m 


2M- 
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and the fraction x' is replaced by the longitudinal momentum k^'^ 


kf'^ = { x' - ^ ) M20 + 


ms — m 


2M 


20 


Explicitly, one has 


Xi(x, k]) = 


^/Af,„[l - (mj - m2)2/M*] 


2(1 ^) \l MIq — (rrii — m)2 


(29) 


(30) 


(31) 


(32) 


with k'^ = k‘j_ + k^. This procedure amounts to a series of reasonable (but naive) guesses about what 
the solution of a relativistic theory involving conhning interactions might look like. 

The matrix element of the vector current Q 2 'y~^Qi between pseudoscalar and vector mesons is 
given by 

=< P2e{+l)\Q2l^Qi\Pi >= j d^k^X2{x',k^l)xi{x,kl)-ll;^\ (33) 

0 


where 


= /ii;i 2 [Tri?]*'+in(p 2 , X 2 h^u{pi, Xi)Poo] 


(34) 
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To define the axial form factors we need two matrix elements of the axial current Q 2 'l^lbQi corre¬ 
sponding to the transverse and longitudinal polarization states £(+l) and £( 0 ) of a vector meson 

r{q'^,q±) 

JA'\q^Q±) =< A^(+1)IQ27+75Qi|Pi >= / v/■ 4+'^ (35) 

0 

=< P2£{0)\Q2'l^'l5Qi\Pi >= J ~ J (^‘^^±X 2 {P, k'l)xi{x, kl) ■ (36) 

0 


where 


Ia^'’ = fiiiJ, 2 Tr[Rt+MP 2 , A2)7’^75m(pi, Ai)i?oo], (37) 

= pi/i2Tr[i?i^oM(p2, A2)7+75 m(pi, Ai)i?oo]- (38) 

Noting that h(p 2 , A 2 ) 7 +m(pi, Ai) = h(p 2 , A 2 ) 7 ’^ 75 m(Pi, Ai) = ^ptP 2 Px 2 ^^P\i, 

where ^px are the Pauli spinors and using the explicit expressions for light-cone spin structure func¬ 
tions -Rjjg from [§, |] we obtain 



/y(a;;r;fc_L;g_L) = [AA + k^kJ, 

( 39 ) 


I^v \x]r]kr,q 2 _) - [k[Ai ^ ], 

x' M20 + rn + m2 

( 40 ) 


I^A \x;r;k±;q^) - [{ 2 x' l)k[Ai + kiA2 + 2 k[ , ^ 

x' M20 + m + m2 

(41) 


, N 2 Mi r ,, ,, , , , (1 — 2a;0M2o + n7.2 — m,, , 

I\\x;r;kr,q±) = , 2x (1 x)M2oAi+^ a# 1 1 (k^k^ + A^) , 

x' M20 + m + m2 

( 42 ) 

with 

Ai = xmi + (1 — x)m, A = x'm2 + (1 — x')m, 

( 43 ) 

and 

B = {1 — x')m — x'm2. 

(44) 


4 The definition of the vector and axial form factors 

We hrst reexamine an approach of ref. to calculate the form factors Fi{q^) for the PS — PS 
transitions. Eq. ( 6 ) for the plus component of the vector current yields only one constraint for the 
two form factors Fi{q^) and Fo(g^). In order to invert Eq. ( 6 ) the matrix element of the current 
was calculated in p] in two reference frames having the 3-axis parallel and anti-parallel to the 3- 
momentum of the daughter meson. This corresponds to the choice ai = 0, 02 = vr where the angle 
a is dehned in Section 2. One might ask what happens when we use two other frames. Specifying 
these frames by the two arbitrary angles ai and 02 we can now write Eq. (35) of ref. 0 as 

(1-r(a2))Jy(g^r(al)) - (1-r(ai))Jy(g^r(a2)) 

> =-mxmTPWT-■ 
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In this equation we can take the limit 02 —> ai = a which would correspond to the use of only one 
reference frame specified by the angle a. In this case the form factor Fi{q^) is expressed in terms of 
the matrix element of the current and its derivative. Indeed, from Eq. (45) one derives 


Fi{q^) 


1 

2ii^ 


Jv'(g^r(a)) + (1 


r(a)) 


dJv{q‘^,r{a)) 

dr 


(46) 


It is understood that when calculating the partial derivative the variable q± entering the integrand 
of Eq. (22) is expressed in terms of and r using Eq. (2). 

The final choice of a can be made by the requirement that at the point of the maximum recoil, 
q^ = 0, Eq. (46) must coincide with the standard LF result for < 0 obtained in the Breit 
frame p, |^. One can easily verify that this condition is fulfilled only for a = 0. Eq. (46) for 
a = 0 coincides with that of refs. [Q, || at g^ = 0 and g^ = g^a^,, but the two definitions differ 
at intermediate g^. Recall |p that the time-like LF result for the “good” component of the weak 
vector current J~^ = coincides with the contribution of the spectator pole of the Feynman 

triangle diagram, which corresponds to the valence quark approximation, while the remaining part 
of the Feynman diagram, the so-called Z graph, can not be expressed directly in terms of a valence 
quark wave function.^ The sum of both contributions does not depend, of course, on the choice of 
the frame but each contribution is frame-dependent. Therefore the time-like LF result for the form 
factors generally depends on the recoil direction of the daughter meson relative to the 3-axis e.g. on 
the choice of angles a specifying a reference frame. Fortunately this dependence is marginal even for 
heavy-to-light transitions. We shall illustrate this point later on. 

Now we apply the same procedure to the calculation of the form factors governing the PS — V 
transitions. The vector form factor R(g^) can be found from Eqs. (7, 33). Note that for the pure 
longitudinal configuration both the integral of Eq. (33) and the covariant 


= -^MiMaV'^^sina = (47) 

vanish. Therefore to calculate the form factor V (g^) we will choose a 7 ^ 0 at the outset and let a —>■ 0 
at the end of calculations to get 


V{q^) 


V2 


(1 + 0 




J qj_=0 


To define the axial form factors we first note that 


h+i) 

•Ja 


1 ,l+r(a) . o. l+ria) , n,, 

-j=Miq^{ -wvj—a+(g ) H--—a_(g )). 

V 2 r[a) r[a) 


(48) 


(49) 


Then we proceed exactly in the same way as was described above for calculating the form factor 
Fi(g^). An additional difficulty is that the integral of Eq. (35) also vanishes for a = 0. Therefore to 
invert Eq.(35) we first solve it for a = ai and a = 02 and then take the limit cxi, 02 —> 0. We obtain 


an7) = ±^W(/lW(r)+r(lq:r)4h^), (50) 

^ Note that the contribution of the Z-graph is expressed in terms of the integral over the region r < x < 1, therefore 
to diminish its contribution one must choose r (at given q^) to have its maximum close to 1; in this way we again 
come to the choice a = 0. 
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where 


yl(+)(r) = 


9 j{+i) 

dq/^ 


JV^\q\q±) 


g±=o 


In order to calculate the remaining form factor f{q^) we use the longitudinal polarization ^(O) 
get 


( 61 ) 


to 


f{q^) = 


(g^ g±) + “ C^)A(+)(r) 


J <jx=0 


At the point of the maximum recoil one obtains from Eqs. (48), (50) and( 52) 


(52) 


1 /( 0 ) 




(1 + 0 


d 

dq± 


4*”(0.9±) 


<?±=o 


( 53 ) 


and 


6lo(0) = ^jf(0,0), 


^ 2 ( 0 ) — 


1 + C 
02 


d 

dq± 


A‘’(0.9x) 


J 5^=0 


( 54 ) 

( 66 ) 


where the relationships (13, 15) between the form factors Aq, A 2 and the form factors a± and / 
have been used. We show in the Appendix that Eqs. (53-55) coincide with the standard result P, |] 
obtained using the Hamiltonian LF formalism in the Breit frameQ. 


5 Results 


The calculation of the semileptonic form factors Fi(g^), V{q‘^), Ao(g^), Ai^q"^) and A 2 {q‘^) has been 
performed using our LF results (Eqs. 46, 48, 50, 52) for computing the matrix elements of the vector 
and axial weak currents. The LF meson wave functions Xi constructed according to the rules 
given in Section 3. For the radial wave functions appearing in Eq. (32), the Gaussian ansatz of the 
ISGW model has been adopted in which the main characteristic is conhnement^ The values of the 
parameters are taken from the updated version [|^ of the ISGW modelQ. We use the physical values 
for the meson masses taken from the latest Particle Data Group publication [p!3| . 

Our results for the form factors at = 0 and semileptonic decay rates T (Eqs. (16, 17)) are 
collected in Tables 1 and 2. The behavior of the form factors Ti(g^) for B —> Dii>£, B —> niui, 
D —>■ Kiui, and D —>■ 7rii>£ semileptonic transitions is shown in Figs. 1-4, respectively, while the g^ 
behavior of the form factors l/(g^), Ao(g^), Ai(g^), and A 2 (g^) for B D*£h'£, B pih'£, D —> K*ih'£, 
and D —>■ pii>£ transitions is shown in Figs. 5-8. Note that near the zero-recoil point y = 1 the 
form factors do not obey at all the pattern of pole dominance and even may have a negative slope. 
This behavior, especially prominent for the b ^ u transitions, can be ascribed to the fact that our 


^We have compared our numerical results for the PS V form factors with those of Jaus |l0l u sing the set of 
quark model parameters (the constituent quark masses and harmonic oscillator lengths) found in [|10| by exploiting 
the duality of the VMD and CQM pictures. This comparison is presented in Table 1. Note that there is a misprint 
in Table VI of ref. [^. The correct Table is obtained if the symbols V, Ai and A 2 are replaced by Ai, A 2 , and V . 
Our results for the vector and axial form factors for B and D decays are very similar to the results of ref. 0 but the 
rates are a bit different and reflect the difference in the dependence of the form factors. 

"‘There is some model dependence here; other choices Q for w{k'^) will lead for somewhat different predictions. 

^In what follows this updated version will be referred as the ISGW2 model. 
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calculations do not include the contribution of the quark-pair creation diagram (the Z-graph) which, 
in another language, corresponds to the contribution of \B*7t > or \D*7r > states [Q. However, 
the corresponding decay rates are not affected by the Z-graph because the contribution to the rates 
arising from the region near y = 1 is kinematically suppressed. 

In order to estimate uncertainties introduced by a particular choice of the reference frame we have 
compared in Table 2 the decay rates for the PS — PS transitions calculated from our new Eq. (46) 
(corresponding to the choice ai,a 2 —> 0) and from Eq. (35) of ref. (corresponding to the choice 

= 0, 02 = vr). We have found the new decay rates to be slightly smaller than those of ref. 0: the 
differences comprise 7% for b ^ c, 9% for b ^ u, 13% for c —> s, and 4% for c —> u transitions. We 
conclude that a good stability of the results is always maintained and the decay rates (with exception 
of c —> s case) are predicted with accuracy better than 10%. The same conclusion is preserved for 
more complicated meson wave functions corresponding to the inclusion of the one-gluon exchange at 
small distances 

We now discuss the results arranged in order of decreasing active quark mass. We will compare 
our results to experimental data and to predictions of different models. 


5.1 Decays B Diug and B D*iiy£ 


First, we want to estimate the values of \Vcb\ based on our predictions of the normalization of the 
form factors and their dependence. From the partial widths^ 


r(H° - 

D-Bu,) = 

(1.22 ±0.32) ■ 

lO'o 

(56) 

r(H° - 

D-*£+i^e) = 

= (2.92 ±0.17) 

■ 101° s-i 

(57) 

T{B^ - 

D-Bue) and T{B° D*- 

Bui) from Table 2, we obtain 


iHcl 

= 0.037 ±0.004, {B^ ^ D 


(58) 

iHcl 

= 0.036 ±0.002, {B^ ^ D- 

-*Bue). 

(59) 


The results are consistent in that the extracted values of Vbc are consistent for two decay modes. 

dq^ 


Fig.9 shows the ^ distribution calculated for |Vbc| = 0.036 and |I4c| = 0.0374. The last value has 


been obtained from our fit to the recent CLEO data 0- 
with the updated “experimental” determinations of iHc 


Both LF predictions are in agreement 


19 


obtained from exclusive and inclusive 
semileptonic decays of the B-meson {\Vbc\exci = 0.0373 ± O.OOdSexp ± 0.0065i/i and \Vbc\inci = 0.0398 ± 
O.OOOSexp ± 0.0040i/i), and also agree within the error bars with the recent determination \Vbc\ = 
0.041 ± 0.003 ± 0.002 from a partial rate corresponding to the kinematical region where 


D* recoil slowly [^. In Table 3 we compare the values of |I4c| for several quark models. 

Our predicted values Ri = I7(0)/Hi(0) = 1.08 and R 2 = H2(0)/Hi(0) = 0.89 are in good 
agreement with the preliminary £t results by CLEO Ri = 1.18 ± 0.30 ± 0.12 and R 2 = 0.71 ± 
0.22 ± 0.07. For the branching ratios BR{B^ BR{B^ —> D*~Bvi) we obtain 


BR{B'' D-rue) = 2.12 




cb 


0.039 


( 60 ) 


®These values have been derived by combining the branching ratios Br{B^ —s- D Pi'i) = (1.9 ± 0.5)% and 
Br{B^ —> D~*Pvi) = (4.56 ± 0.27)% with the world average of the B^ lifetime t^o = 1.56 ± 0.06 ps |^|. Our 


estimation of r(i?° ^ agrees with the recent result by CLEO collaboration r(i3° ^ D~*Pvg) = 

[2.99 ± 0.19{stat) ± 0.27{syst) ± 0.20{lifetime)] ■ 10^° based on the average of the LEP and CDE measurements 
of the lifetimes j^, rgo = 1.53 ± 0.09 ps, tb+ = 1.68 ± 0.12 ps. 
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BR{B^ D*-tui) = 5.40 




0.039 


(61) 


These values only slightly exceed the corresponding branching fractions obtained using the heavy- 
quark symmetry approximation i.e. adopting a single Isgur-Wise function for all the form factors: 

BR{B^ —>• = 1.65 , and BR{B^ —>• = 5.17 . Our prediction for the 


, 0.039 , 

ratio of of the longitudinal and transversal contributions to the B 

Tr 


with the recent result by CLEO Collaboration ^ = 1.24 ± 0.16 


r{B^ D*-i+Ui) 
r (50 ^ D-i+Ui) 


= 2.54, 


D* decay is 1.17 which agrees 
We also obtain the ratio 

(62) 


that agrees with experimental result 2.4 ± 0.6 |]T^. Some other model predictions for the B ^ D* 
decay are collected in Table 4. 


5.2 Decays B Trivi and B piui 

We now consider the semileptonic B-meson decay corresponding to the quark process b uiui. The 
investigation of this decay is very important for the determination of the CKM matrix element 114.^1, 
which plays an important role for the CP violation in the Standard Model. The determination of Vbu 
is one of the most challenging measurements in B physics. Very recently, the CLEO Collaboration 
M nas reported the first signal for exclusive semileptonic decays of the B meson into charmless 
hnal states, in particular for the decay modes B —> Tiiui, B —> However, there is a signihcant 

model dependence in the simulation of the reconstruction efficiencies. The observed branching ratios, 
assuming the efficiencies obtained from the ISGW model are 


Br{B TTiui) = (1.34 ± 0.45) ■ 10-^ Br{B piue) = (2.28l°|^) ■ 10 


-4 


(63) 


From our predicted rates r(i?° 
(see Table 2) we obtain 
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bu 


Vcb 

Vbu 


V 


cb 


= 0.087 ±0.015, {B 
= 0.092 ±0.015, {B 


l2 — 

ps 

^ and r(H° - 

p ^ 13.2 ps ^ 

(H° 

1 

+ 

(64) 

(H° 

it 

^ p Bua). 

(65) 

(H° - 

7r"£+z/£), 

(66) 

(H° - 

p~Bue)- 

(67) 


These values are in agreement^ with the value derived from measurements of the end-point region of 
the lepton spectrum in inclusive semileptonic decays ^ 


VIZ. 


Vbu 


Vcb 


incl 


= 0.08 ± O.Olexp ± 0.02t/i. 


( 68 ) 


^If the BSW model ^ is used to reconstruct efficiencies the branching fractions are somewhat higher: 

Br{B —> ttRi) = (1.63 ± 0.57) • 10“"^, and Br{B —> pivi) = (S.SSt^'gg) • 10“^. Combining the average of experimental 


results for the ISGW and BSW models with our predicted rates one obtains |T4„| = 0.0033 ± 0.0004, 
0.015, (B° ^ TT-f+i/^ and |14« I = 0.0039 ±0.0004, ^ = 0.107 ± 0.015, {B° ^ p-^+vi). 




Vcl 


= 0.092 ± 
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A large variety of calculations of the ratio | Vb^i/Vfed, based on various non-perturbative approaches, 
exists in the literature; the results typically lie in the range from 0.06 to 0.11. A sample of these 
calculations is given in Table 5. 

Note that our result for the ratio of branching fractions 


r(i;° ^ p-e*n) 

T{B^ —>• TT~i+iy£) 


compares favorably to the recent CLEO result 


0 


T{B^ ^ p-£+z/,) 
r(S 0 ^ TT-i+ue) 


1.70d°|° ± 0.58 


+0.00 

-0.34 


(69) 


(70) 


The dispersion relation approach to the ISGW2 model yields for this ratio the values of 1.17—1.34. 
Our prediction for the ratio Tl/Tt for the B —> piui decay rate is 0.92; other models yield Tl/Tt = 


1.13 


1.34 


0.3 


], 0.82 ig, and 0.52 ± 0.08 


5.3 Decays D Kii/£ and D K*ih'£ 


In the charm sector, the absolute scale of theoretical predictions for the form factors which char¬ 
acterize the hadronic D K* and D ^ p matrix elements of the weak currents can be tested, 
because the CKM elements can be determined independently of the D semileptonic decay rate us¬ 
ing the assumption of CKM unitarity and the smallness of the CKM matrix elements for B decay 
^]. Semileptonic decays D —>■ £ = e, p have been studied extensively and the form factors 


have been determined (with some assumptions concerning their shape). In particular, the CLEO 

K{K*)£i/i decays in one experiment: 


collaboration ||^ has measured all four D 

T{D K£ui) = (9.1 ± 0.3 ± 0.6) ■ lO^V^ 
T{D K*lvi) = (5.7 ± 0.7) ■ 10^V\ 


R = 


T{D ^ K*lvi) 
T{D ^ Kiue) 


Using \Vcs\ = 0.974 we predict 


= 0.62 ±0.08. 


T{D Kiui) = 9.38 ■ 10^°s-\ 

T{D = 6.36 ■ 10^°s-\ 

7 ? = 0 . 68 , 

in an agreement with experimental data. Some other model predictions for the D 
collected in Table 6. 


(71) 

(72) 

(73) 

(74) 

(76) 

(76) 

K* decay are 


5.4 Decays D Trii/i and D pii/i 


Assuming 114^1 = 0.221, which is inferred from the unitarity of the CKM matrix 
the semileptonic decay rates the following values: 


131, we predict for 


r(T)° ^ vr-e+i/e) = 7.4 ■ lO'^ ps-\ 


(77) 
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T{D^ —>• p e^Ve) = 3.8 -10 ^ ps 

These predictions should be compared with the experimental result 

^ vr-e+i/e) = (9.4l|5) • ps-\ 


and the recent light-front QCD sum rule prediction 


Tsr{D<^ ^ vr-e+i/e) = (7.6 ± 0.2) ■ ps-\ 


The CLEO collaboration recently determined the following ratio of decay rates 

r(T)° ^ TT-e+Z/e) 


r(T)0 ^ K-e+Ue) 

while the 77653 collaboration [^| measured 

r(D 0 ^ 

r(T )0 ^ K*-e+iye) 


= 0.103 ±0.039 ±0.013, 


(78) 

(79) 

(80) 

(81) 


= 0.088 ±0.062 ±0.028. 


(82) 


Our calculated fractions are 0.079 and 0.060, respectively, which seem to be a bit small; nevertheless 
they are in agreement with the experimental values within the large error bars. The dispersion 
approach to the ISGW2 model pT[ yields for these fractions 0.071 and 0.047 — 0.056. 


6 Conclusions 

In this paper we have examined the weak transition form factors which govern the heavy-to-heavy 
and heavy-to-light semileptonic decays of pseudoscalar mesons within a relativistic constituent quark 
model based on the light-front formalism. We have obtained all the relevant expressions for the weak 
vector and axial form factors. For numerical estimates we employed the LF wave function, that is 
related to the equal-time wave function of the ISGW model. The transition form factors have been 
calculated in the whole kinematical region accessible in semileptonic decays. We have calculated 
the form factors and the decay rates for the B D[D*)iPe, B 7r(p)f'z/^, D —^ K{K*)£i>£ and 
D —> n{p)£pi semileptonic decays. Our results have been successfully compared with available 
experimental data. In particular for b ^ c and b ^ u transitions, using the available experimental 
information on the semileptonic decay rates, the GKM parameters have been estimated. 

Since our calculations neglect the contribution of the pair creation from the vacuum it would 
appear that the good agreement with the data suggests that these might very well be small. However, 
an estimate of such contributions, particularly in case of the heavy-to-light transitions, is important 
for there to be a complete comparison with experimental data. 


Acknowledgment 

We thank I.L. Grach for the collaboration in the early stages of this work. The authors acknowledge 
the financial support of the INTAS-RFBR grant, ref. No 95-1300. This work was in part supported 
by the RFBR grant, ref. No 95-02-04808a and by the Natural Sciences and Engineering Gouncil of 
Ganada. 


12 













Appendix. Form factors at = 0 


In this Appendix we show that Eqs. (53-55) for V (0), Ao(0) and ^ 2 ( 0 ) coincide with the ones obtained 
using the LF formalism in the Breit frame for < 0. Using our notation we write^ Eqs.(4.4, 4.5, 
4.6a, 4.9) of Ref. ^ in the form 

^ I ^ I A:±;gx), (A.l) 

^ ^ 0 

^ 2W I ^ I ^^^^^ 2 ( 3 :, A;(^)xi(a;,^l)2^A^(a;;^±;?±), (A.2) 

^ 0 

^ I ^ I ^^^^^2(a:,/c2)xi(a^,fci)2i^^^(a^;fc±;g±), (A.3) 


where 


Iv^\x]kj_]q±) = y/2Mi 
I^A\x-,k±;q±) = 2M, 


mi-m 2 , 

Ai - ^ —k^q± + 


q^ M20 + m 2 + m 

(1 — 2x)M2o + m2 — m 


k] + 




2(1 x)M2oAi+ 

{x; k±; q±) = —V2Mi [{2x — l)^i- 

, X _ Ekai 

"‘"^"*' -[2(1 — x)m + m 2 + {2x — l)mi] + 2- — - 


(k^k'^ + AiB) 


Q 


M 20 + m 2 + m 


(k^kl + A2B) 


(A.4) 

(A.5) 

(A.6) 


and k^ = k_i_ -|- xqj_ . 

Noting that in our kinematics the momentum transfer is anti-parallel to the 1-axis, i.e., q = 
(g°, q^, g^, g^) = (0, —q±, 0, 0) we obtain after simple algebra 


Iv^\x] k±] q±) = —ll^^\x, 1, k±, q±), 

(A.7) 

Xf^(a;; k^; g^) = 4°^(x; 1; k±; g^). 

(A.8) 

Ia^\x; k±; qx) = —Ia^\x; 1; k±; qx), 

(A.9) 


q± 

In the limit g± —> 0 using Eqs. (33, 35, 36) one easily derives that Eqs. (A.l, A.2, A.3) coincide with 
Eqs. (53, 54, 55). 


®Note that our variables and notation are related to those of ref. Q by 

p([ = — k'j_, p'^ = — k^, ^ = 1 — a;, m 2 = m, m" = m 2 , m( = mi, 
The quantity n(p",p') defined by Eq.(4.2) of ||] can be written in terms of our LF functions as 

^-^—n{kkk) = -X2{x,k'l)xi{x,k'i) 

ax I — X X 
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Table 1. Form factors at = 0 and decay rates for the set of quark model parameters of . 
The decay rates for be and bu transitions are in the units \Vbc\‘^ps~^ and \Vbu\^ps~^, respectively. 
The decay rates for cs and cd transitions are in units The values of the CKM parameters 

\Vcd\ = 0.22, |14s| = 0.9743 are used. In parentheses are given the numerical results of ref. 


Q 1 Q 2 

n(o) 

1/(0) 

Ai(0) 

^ 2 ( 0 ) 

r(ps PS) 

r(PS V) 

be 

0.68(0.69) 

0.79(0.81) 

0.67(0.69) 

0.61(0.64) 

9.3(9.6) 

24.5(25.3) 

bu 

0.25(0.27) 

0.34(0.35) 

0.26(0.26) 

0.24(0.24) 

8.1(10) 

18.3(19.1) 

cs 

0.78(0.78) 

1.03(1.04) 

0.65(0.66) 

0.43(0.43) 

9.7(9.6) 

6.2(5.5) 

cd 

0.66(0.67) 

0.90(0.93) 

0.57(0.58) 

0.42(0.42) 

0.7(0.8) 

0.39(0.33) 


Table 2. Form factors at = 0 and decay rates for the set of quark model parameters of . 
The decay rates for be and bu transitions are in the units \Vbc\'^ps~^ and \Vbu\'^ps~^, respectively. 
The decay rates for es and eu transitions are in |l/cspl0^‘’s“^ and respectively. In 

parentheses are given the decay rates for the PS — PS transitions taken from ref. P] . 


Q1Q2 

n(o) 

1 /(0) 

Ai(0) 

^ 2 ( 0 ) 

/ 4 o( 0 ) 

r(PS PS) 

r(PS V) 

be 

0.683 

0.677 

0.623 

0.556 

0.678 

9.09(9.78) 

23.10 

bu 

0.293 

0.216 

0.170 

0.155 

0.214 

8.72(9.62) 

13.2 

cs 

0.780 

0.777 

0.633 

0.464 

0.73 

9.88(11.30) 

6.70 

cd 

0.681 

0.663 

0.502 

0.366 

0.60 

15.29(16.0) 

7.85 


Table 3. Model-dependent predictions of the B —> D*ijy partial width Va values derived from 
the measured partial width. An additional 3.5% systematic error due to B lifetime measurements is 
common to all the \Vcb\ values given. 



This work 

ISGWg 

BSWg 

Ks|g 

r(B^n*£i2)IKbl-'^ 

23.1ps-^ 

24.6ps“^ 

21.9ps“^ 

25.8ps-^ 

lo^-IKbl 

35.6 ±1.1 ±1.6 

34.8 ±1.1 ±1.6 

37.5 ± 1.2 ± 1.7 

34.4 ± 1.1 ± 1.5 


Table 4. Form factors at g^ = 0 and the decay widths in for the semileptonic B^ D* transition. 
The model predictions for the width are in iHcP ps~^. 


B - 

D,B>* 

V'(0) 

i/(0)Mi(0) 

A2(0)/Ai(0) 

r(B*) 

F(T))/F(D*) 

Tl/Ft 

This work 

0.68 

1.08 

0.89 

23.1 

0.38 

1.15 

DR 

21 



0.68 

1.08 

0.89 

21.0 

0.41 

1.28 

ISGW2P1 




24.8 

0.48 

1.04 

Ksen 



1.54 

1.39 

25.8 

0.32 


BSW 

23 


0.71 

1.09 

1.06 

22.0 

0.36 


Stech 

32 


0.75 

1.10 

1.02 

24 

0.37 


SR li 



0.58 

1.26 ±0.08 

1.15 ±0.20 

17 ±6 

0.53 ±0.11 


Exp. 

18 



1.18 ±0.30 ±0.12 

0.71 ±0.22 ±0.07 


0.46 ±0.2 

1.24 ±0.16 
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Table 5. Values for |^| extracted from CLEO measurement of exclusive semileptonic B decays 
into charmless final states. Following Neubert we take in this comparison |14fe| = 0.040. An 
average over the experimental results obtained using the ISGW and BSW models to reconstruct 
efficiencies is used for all except the ISGW and BSW models, where the numbers corresponding to 
these models are used. For our LF model we use the values quoted in Eq.(64, 65). The hrst error 
quoted is experimental, the second (when available) is theoretical. For the recent LG SR result 
only the theoretical error is quoted. 


31 


Method 

Reference 

B —>■ Tliu 

B piu 

LF QM 

this work 

0.078 ±0.013 

o.o837tJ:(;i^ 

QGD sum rules 

Narisc 

Ball 

KhR 

)n 

31 

M 

35 


0.159 ±0.019 ±0.001 

0.085 ±0.010 

0.066^^;^^^ ± 0.003 

0.086 ±0.013 

Lattice QGD 

UKQGD0 
APE H 

0.103 ±0.012+^;^{^ 
0.084 ±0.010 ±0.021 


pQGD 

Li and Yu 

^9 

0.054 ±0.006 


Quark models 

BSW 
KS d 
ISGW 

pa 

3 

1 

11 


0.093 ±0.016 

0.088 ±0.011 

0.074 ±0.012 

0.0767jj:(;i^ 

0.056t°:“® 

o.o79«:°12 


Table 6. Form factors at = 0 and the decay widths in 10^° s ^ for the semileptonic K* 

transition. For our calculations we use |14s| = 0.974 


K* 

V(0) 

Ai(0) 

^2(0) 

r(A:Q 

rL/rr 

This W( 
DR 21 
SR p6 
Lat. 3 
Lat. 3 
Stech 
Exp. jl 

3rk 

8 

7 

321 

3 

0.777 

0.777 

1.1 ±0.25 

1.08 ±0.22 

1 Dl 

1.07 

1.1 ±0.2 

0.633 

0.633 

0.50 ±0.15 

0.67 ±0.11 

n yn+o.or 

U- ' U_o.io 

0.69 

0.56 ±0.04 

0.464 

0.464 

0.60 ±0.15 

0.49 ±0.34 

0 . 60 l°:J° 

0.73 

0.40 ±0.08 

6.36 

5.38 

3.8 ±1.5 

6.9 ±1.8 
6 .0l?| 

7.1 

5.7 ±0.7 

1.30 

1.31 

0.86 ±0.06 
1.2 ±0.3 

1.06 ±0.16 

0.97 

1.15 ±0.17 
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Figure Captions 

Figures 1, 2, 3 and 4 The form factors Fi{q^) for and D ^ tt transitions, 

respectively, obtained nsing the ISGW parameters, bnt adopting the experimental valnes for 
the meson masses. The relation between the kinematical variables y and is given by Eq.(4). 
The solid lines are the results of our LF calculations obtained using Eq.(46) with ai,a 2 —^ 0. 
The dashed lines are the previous results of ref. [Q corresponding to the choice of the two 
reference frames with ai = 0, a 2 = as explained in the text. The corresponding decay rates 
are listed in Table 2. 


Figures 5, 6, 7 and 8 The vector V{q^) and axial form factors AQ{q^), Ai{q^) and A 2 {q^) for 
B —> D*, B ^ p, D K* and D ^ p transitions, respectively, obtained using the ISGW 
parameters, but adopting the experimental values for the meson masses. The solid, dashed, 
dotted and dot-dashed lines represent the vector form factor V{q‘^) and axial form factors 
Ao(g^), Ai{q^) and A 2 {q‘^), respectively. The corresponding decay rates are listed in Table 2. 


Figure 9 The ^ distribution for B —> D*ih'g decays. 


corresponds to our £t to the GLEO data 
0.036 found from the partial width. 


161 while 


The upper thin line with |14c| = 
the lower solid line corresponds to 


0.0374 

iHcl = 
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Figure 2: Form factor Fi{q^) for B ^ tt 
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Figure 9: The distribution for B D*iui decays. 
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